Magneto-transport and divergent screening of driven two dimensional electron gas 
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In two dimensional electron system in magnetic fields such that the Fermi energy lies in high 
Landau levels and driven by microwave ac-field a divergence of the Coulomb screening due to the 
Ryzhii photocurrent is predicted. A model of magneto-transport with superposition of long range 
and short range disorders is introduced. In this model the larger is the gradient of the long range 
potential the smaller is the longitudinal conductivity /resistivity. As a consequence a qualitative 
theory of the recently discovered zero-resistance states is given. 

PACS numbers: 72.10.-d, 73.40.-c, 78.67.-n, 73.43.-f 
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I. INTRODUCTION 

Electron systems interacting with microwaves could be 
driven into highly non-equilibrium state with unusual 
properties. For example two dimensional electron sys- 
tems (2DEG) in GaAs heterostructures subjected to a 
weak perpendicular magnetic field show microwave in- 
duced resistance oscillations which in the limit of ex- 
treme power are capped by zero resistance states^. 
Ryzhii has predicted that an additional electric current 
will flow against or along applied electric field in the 
2DEG driven by microwaves depending on the ratio of 
microwave frequency u> and the cyclotron frequency ujh- 3 
Recent works^ have clarified the origin of the Ryzhii 
current and resistance oscillations in disordered systems. 
At strong microwave driving the negative Ryzhii cur- 
rent could overcome usual positive dissipative current 
and thus results in the total local negative conductiv- 
ity. Negative conductivity leads to the current domain 
instability and the domains could arrange themselves to 
give a zero total resistance^. 

Zero resistance states may have different origins. For 
slow driving u> <C ll>h a wide suppression of resistance^ 
could be a manifestation of 'driven' degeneracy of elec- 
tron states^. Fast driving u> ^> luh averages out the dis- 
order potential and would lead to an ideal state with zero 
resistance (not reported yet to my knowledge) , similar to 
the induced transparency in optics. This paper is con- 
cerned with the intermediate frequencies lj ~ luhi where 
multiple zero resistance states specific only to 2DEG in 
magnetic field emerge due to the Ryzhii phenomenon^. 

In this paper I study effects of the long range disorder 
due to the Si-doped layer in high mobility heterostruc- 
tures. If a long range potential is present in the elec- 
tron system driven by microwaves then random Ryzhii 
currents will be locally induced even in the absence of 
the probe cic-voltage. Electrons will move away from 
the equilibrium density distribution, thus, modifying the 
Coulomb screening. The relationship between the bare 
long range potential Uo(r) and the screened potential u(r) 
is linear one: u(q) — F(q)uo(q), where F(q) is the screen- 
ing factor, provided both \u\, \uq\ <C ep. One of the cen- 



tral results of this paper is that: 



<jd{W)J 



where <jd{W) is a diffusion part of the total conductivity 
a(W) = a D (W) — <7r(W) and a R (W) is the Ryzhii con- 
ductivity (assuming that ku>H < w < {k + l/2)uJH, where 
k is an integer). They both depend on the intensity of mi- 
crowaves W. In the long range limit qR c <C 1, where R c 
is the cyclotron radius, the electron transitions amount 
to a continuous current which is the linear response to 
two different perturbations: weak external electric field 
and weak gradient of electron density: 



ej(f) = a{W)Vu{r) + e 2 D(W)Vn{r) 1 



(2) 



where D(W) is the diffusion constant. The Einstein rela- 
tionship states the detail balance at equilibrium: er(0) = 
e 2 J/(ei?)£>(0) (ctr(O) = 0). It is violated in the driven 
state with cr D (W) = e 2 v(e F )D{W). The source of elec- 
trostatic potential u(r) (in this paper it means the elec- 
trostatic energy) is either donors or electrons: 



u(q)=u {q) + (2ire 2 /\q\) n{q). 



(3) 



The total current under the stationary microwave power 
is zero j(q) = 0. Assuming that electron relaxation due 
to phonons induces no additional current we find from 
Eas. H2l3(l the screening factor Eq.JI}. In this paper a 
connection between the divergence of the screening factor 
Eq.JTJ at some critical power W c and the development 
of a zero resistance state is sought. Toward this end a 
realistic model of magneto-transport without (and then 
with) microwaves in GaAs heterostructures is required. 

There are two well understood models of the magneto- 
transport of 2DEG in high Landau levels. First is the 
model with short range disorder and large scattering mo- 
menta: qln 3> 1, where Ih is the magnetic length - the 
self-consistent Born approximation, exact in high Landau 
levels^. It predicts the classical Drude conductivity ten- 
sor specified by the scattering time r in weak magnetic 
fields loht *C 1 and oscillating longitudinal conductivity 
with: (T™ ax (i/) - e 2 N/h^l where TV is the Landau level 
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index, in strong magnetic fields lort 1. Second is the 
model with long range disorder and small scattering mo- 
menta: qR c <C l,- 11 It employs the approximate conserva- 
tion of the adiabatic invariant and predicts the classical 
localization 12 with the negative magneto-conductivity: 
&xx{H) ~ exp(—H/H c )m& In GaAs heterostructures the 
disorder potential from Si donors has a cut-off in momen- 
tum space q c ~ l/2d where d is the spacer width. Thus, 
a crossover from the short range disorder at weaker mag- 
netic fields to the long range disorder at stronger mag- 
netic fields occurs at Ih ~ 2d or H c ~ 0.3kG (for typical 
d = 80nm). But above this H c experiments^ show pos- 
itive rather than negative magneto-resistance. Second 
contradiction is the behaviour of the Shubnikov-de-Haas 
oscillations when the longitudinal resistivity approaches 
zero. The theory of long range disorder predicts the clas- 
sical localization and wide platoes of quantum Hall zero 
resistivity whereas experiments 2 *^ show no platoes but 
instead a cusp-like touch of zero resistivity value. 

In this paper a model of magneto-transport in high 
mobility GaAs heterostructures with disorder that is a 
superposition of weak short range and strong long range 
potentials is introduced. The main result is that the lo- 
cal diffusion constant is inversely proportional to the lo- 
cal electric field E of the long range potential. Indeed, 
this field lifts the degeneracy of the Landau level and im- 
part a velocity and a linear dispersion to the electron: 
e(p) = el 2 H E x p. Therefore the density of states is in- 
verse to the electric field: v{tp) ~ l/E (neglecting a 
logarithmic factor). The mean scattering rate 1/r is pro- 
portional to the mean square of the short-range potential 
w. 1/r = 2'Kv(tF){w 2 ) , in analogy with the usual disor- 
dered metals. In magnetic field the diffusion constant is 
proportional to 1/r due to dominant Lorents force. 

In this paper the temperature is zero as it is appro- 
priate for physics of disorder. It is assumed that in the 
state driven by microwaves the electron relaxation due to 
phonons is effective and keeps the electron distribution 
stationary, approximated by the Fermi-Dirac function. 

In section II a model of free electrons moving in disor- 
der potential is introduced. In section III the diagram- 
matic method is described and the density of states and 
the conductivity tensor for three simplified disorder mod- 
els are calculated. In section IV the magneto-transport 
of 2DEG in the realistic model with both long and short 
ranged disorder is found. In Section V the state driven 
by microwaves is considered and the divergent screen- 
ing Eq.Q is derived in the microscopic theory. It is ex- 
plained how zero resistance states emerge. 



II. MODEL 

In the Fermi liquid theory the transport is described 
in the approximation of non-interacting quasiparticles. 
In strong magnetic fields and for long range poten- 
tial this pictures requires explanation. In 2DEG with 
weakly broaden Landau levels in the long range disor- 



der the chemical potential could not be a smooth func- 
tion. Instead the Coulomb interaction makes the den- 
sity a smooth function whereas the chemical potential is 
pinned by the Landau level and makes jumps along lines 
where the density crosses an integer filling values. It has 
been proven that the width of the quantum Hall stripes 
of integer filling is negligible in high Landau levels^* (for 
H = OAkG and n = 3 x 10 n cm~ 2 the stripe width is 
a = 0.3d <C 2d). In a superposition of long range and 
weak short range potentials a quasiparticle will diffuse 
inside a band around the energy level contour of the long 
range potential 1 -. The boundaries of this diffusion band 
are determined by the conservation of the total energy 
that is divided between the long range potential and a 
'reservoir' of the short range energy (local Landau level 
broadening)^. But interacting quasiparticles diffuse out 
of these band limits. One quasiparticle could climb up 
in the long range potential landscape whereas another 
quasiparticle falls down - very similar to the backflow in 
Fermi liquids. Therefore the magneto-transport in high 
Landau levels derives from an unrestricted quasiparticle 
diffusion. It destroys the classical localization and could 
be described in the non-interacting approximation. 

The Hamiltonian of free electrons in a random disorder 
potential u(r) reads: 



2m \ 



u(r) 



ip{f)d 2 



(4) 

where m is the renormalized by interaction band mass, 
A(r) is the vector potential of the perpendicular mag- 
netic field H. We use the Landau gauge A x — —Hy and 
suppress the spin index due to the vanishing exchange 
interaction in high Landau levels. 

There are three physical parameters of 2DEG: m, H 
and the density of electrons: n = (ip^(r)tp(r)) = N/2nl 2 H . 
One is the Fermi energy, whereas the two others are di- 
mensionless: the Landau level number N and the interac- 
tion parameter: r s = e m/s/irn. In spirit of Landau the- 
ory of Fermi liquids the model J3J describes quasiparticles 
near the Fermi level. We use the magnetic units: h = 1, 
e = c, H = 1, uh = l/m. In high Landau levels (N ^> 1) 
there a re t hree distinct lengths: the cyclotron radius 
R c = \/2NIh, the magnetic length Ijj = ^Jch/eH = 1 
and the inter-electron distance 1/y/im. 

A random electrostatic potential in the quantum well 
is created by remote ionized Si donors with the positive 
charge e. Let donors be confined to a narrow 2D layer 
- (5-doped layer - parallel to the quantum well and sepa- 
rated by the spacer of width d free from impurities. The 
positions of donors assumed to be uncorrelated and their 
areal density is equal to the density of electrons nf® p = n. 
The correlation function of the donor electrostatic poten- 
tial uo(f) reads: 



S (q) = (u (q)u (-q)) 



{2Tre 2 



V \Q\ 



nexp(-2\q\d), (5) 



in the momentum range q <C \fn. We assume dy/n ^> 
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The donor random potential u(f) is Gaussian with the 
probability distribution functional: 



V[u] = iexp(-5>($H-$)/2S(g)), 



(6) 



where Z is the normalization factor. 

Electrons screen the donor potential uo(r) near the 
2DEG plane according to Eq.(j2|. In high Landau lev- 
els a linear response occurs even in the quantum limit of 
non-overlapping Landau levels. The correlation function 
of electron density in response to uo(r) is: (n(q)n(—q)) = 
n exp (~2\q]d) . The average density variation: ((Sn) 2 ) = 
n/2ir{2d) 2 . For typical densities n in GaAs heterostruc- 
tures the screening is complete: 

S(q) = (u(<?)u(-q)) = nexp(-2\q\d)/v 2 (e F ). (7) 

This result is quasi-classical in nature but holds for high 
locally non-overlapping Landau levels as well. 

Usually the areal density of Si in <5-doped layer is much 
larger than n 2 ^ p . Then ionized donors can adjust posi- 
tions to lower the total electrostatic energy. Resulting 
configuration of ionized donors is different from the un- 
corrected one and is better described by highly corre- 
lated liquid state. The correlation function is modified: 



Stf) = nS f (q)exp(-2\q]d)/v 2 (e F ), 



(8) 



where Sf(q) is the structural function of the liquid state. 
Generally, Sf(q) = Aq 2 at small q. We note that for 
periodic structures of ionized donors S(q) — for all q 
except on the Brave lattice points. 

We expand the electron Green function into pertur- 
bative series of the disorder potential u(f), using a dia- 
grammatic method^. We then average this series over 
the probability distribution density Eq.©. In the result 
we get diagrams with electron and impurity lines. The 
electron density vertex describes an electron scattering 
off the potential u(q) from the state (n,p) into a state 
(n',p'), where the indices p and p' counts degeneracy in- 
side the Landau level n. The vertex is the product of the 
reduced vertex and the magnetic phase factor: 

V(np,n'p',q) = V nn > (q)U{p, q) V,p+% • ( 9 ) 

The magnetic phase factor is universal for all Landau 
level: U(p, q) — exp(iq x (p + q y /2)), whereas the reduced 
vertex depends on the incoming and outgoing Landau 
level indices and the scattering momentum: 



V„ 



■ ± iqz 

V2 



-\n — n | 

J min(nn') 



2\ 



(10) 

where L^ l (x) is the Laguerre polynomial, and where the 
sign ± in corresponds to the cases: n > n! and 
n < n' . In the limit n,n' — > oo the quasi-classical ap- 
proximation applies provided q < Ay/nn: 



V nrl '(q) = 



q x + iq y 
\q\ 



qpnn'(q) 



J\n-n'\ (4>nn'{q)) : 
(11) 



where p nn > (q) = + 
momentum and i 



- 1 — f? 2 /4 is the quasi-classical 

(q) = qpnn>(q)/2 + (n + ri + 1) 

arcsin(g/2v^ + n' + 1) is the quasi-classical phase. 

The electron propagation in static potential is either 
retarded or advanced. The bare Green function is diag- 
onal in the Landau level index: 



R.A 



(e) = l/(e - nu H ± iS), 



(12) 



where < 6 — > is the bare dissipation. In the Schwinger 
proper time representation: 



HZ® 



exp (i(e — nu>H)t — 8t) dt. 



(13) 



The advanced Green function is the complex conjugate 
of the retarded one: <?^(e) = (p^(e)J . 

We include two electron density vertices into the im- 
purity line that connects them. Thus the impurity line 
that is not crossed by other impurity lines reads: 



u k (N) = / S{q)\V N , N+k(q)f 



d 2 q 



(14) 



In this notation the Green functions and the impurity 
lines does not depend on the degeneracy index p. 

If the two impurity lines with wavevectors q\ and q\ 
intersect each other then the total integrand given by the 
product of the two integrands in Eq. (|14fl is entangled by 
a complex magnetic phase factor: exp(iq\ x q 2 ). This 
magnetic phase factor can be conveniently represented 
as a commutator: 



exp (iq\ xq 2 )= V M {qi) Vm(£;) 
of the two magnetic phase operators: 

Vfcf(g) = exp (i(qa) + q*a)/V2 



(15) 



(16) 



that substitutes the magnetic phase in Eq.© and is ex- 
pressed in terms of the Bose operators a* and a. This 
operator acts in the space of the Landau level degener- 
acy index p in the Landau gauge. 



III. DENSITY OF STATES AND 
CONDUCTIVITY 

In this section we consider three artificial models of 
disorder using two cutoff wavevectors: qh *C 1/Ih and 
qs ^ l/iff- A short range potential model S is the 
model (0J with the correlation function S(q) — for 
q < qs- A long range potential model L is the model 
Q with the correlation function S(q) = for q > q^. A 
critical model C is the model @ with the special correla- 
tion function: S(q) = A/q, for all wavevectors from 1/R C 
to yfwn. In high Landau levels around TV the correla- 
tion function S(q) is represented by the set of coefficients 
Uk(N) Ea. (|14f> . We expand the average retarded Green 
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matrix G R n ,(e) into a series of coefficients Uk, using a di- 
agrammatic methodic. The density of states is related 
to the imaginary part of this Green function 16 : 



-ImTr(G fl (e))/7r. 



(17) 



Density of states for non-overlapping Landau levels 
in the limit ^/uq/loh -C 1 and Uk <C uq for k > 0. Ac- 
cordingly we set Uk = for k > 1 and uq = u. Thus the 
Green function is diagonal in the Landau level index N 
and we omit it from the general expansion: 



(G R (e))=g R (e)J2c k u k (g R (e)) 



21,- 



(18) 



k=0 



where the combinatorial coefficients Cfc, with Co = 1, are 
determined by the magnetic phase factors and g R (e) is 
the bare Green function (|12[1 . In the limit N — ► oo there 
are just three such unique series c k corresponding to the 
three different models: C, L, S. 

In model S the magnetic phase factor for diagrams with 
crossed impurity lines is small as 1/N. The self-consistent 
Born approximation neglects crossing and coefficients c k 
are the Catalan numbers (1,2,5,14,42...) and the aver- 
aged Green function reads& 



(G R (e))s = 



2ug(e) 



(19) 



where <? _1 (e) = e — Nu>h- 

In model L the the magnetic phase for crossing of two 
impurities lines is zero. Therefore the coefficient Ck is 
given by the combinatorial number of all pairing of 2k 
points: Cfe = (2k— 1)!!. The average Green function 118|) 
has been summed up in the Ref^£: 

r°° ( ut 2 t \ 

(G R (e)) L = i ^ exp (-— +z— J dt. (20) 

For critical model C each integral over impurity line 
wavevector gives a logarithmic factor log(i? c v / 7rn). In the 
leading logarithm approximation we have calculated the 
nine first coefficients Ck and they coincide exactly with 
that of the Wegner model for 'delta'-correlated impurity 
potential on the lowest Landau level 18 : 



(G R (e))c 



d_ 
d~e 



log 



exp 



ut 2 



•i-prjdt, (21) 



despite the fact that the model C is defined in high Lan- 
dau levels N > 1. 

The density of states Ea. (|17(l is the semi-circle law 
in the case of the model S, the Gaussian nonuniform 
broadening in the case of the model L^£, and the Wegner 
distribution 18 for the critical model C. 

Arbitrary Landau level mixing. For long range 
model L we prove Ea. l|2U[l in the case of arbitrary Lan- 
dau level mixing. This subsection essentially rederives 
results of Refiii using an alternative method. For model 



L the electron Green function can be represented in the 
Schwingcr proper time as: 



ut \ - 

l€t — y ~ ^ 

G'L. = 1 1 dt U nn ,(t) e ! 



' 2 d^h 

2tt ' 
(22) 

where u — uq, k > 1 and the evolution matrix: 



U rm i(t) = (n\exp(-iht)\n'), (23) 
is the unitary evolution operator of the Hamiltonian 

h = -iu JH —+Y J ^{^ + r k e-^), (24) 
^ fc=i 

where are complex parameters. This Hamiltonian 
accounts for the Landau level mixing at arbitrary val- 
ues of Uk/uJn- The first term in l|24|l is the kinetic en- 
ergy whereas the non-diagonal second term accounts for 
scattering transfer between Landau levels. The conver- 
gent integrals over complex variables ipk (IV'fcl ~ 1) gi ye 
a combinatorial number of crossed line pairing and the 
term ut 2 /2 generates scattering processes that conserve 
the Landau level index. The wave function 



I / 

Xj(4>) = % ex P I 2J 



ik 



k(- 



%4> 



-i<f>\ 



(25) 

where Z is the normalization factor is the solution of the 
Schroedinger equation hx(<j>) — e x( ( t ) ) with eigenvalue 
e = Jluh- The condition that the wave function is single 
valued restricts j to be an integer. 

We prove that eigenvalues of Hamiltonian l|24l) are jujji 
exactly. Consider shift operators: §k = exp(ik<p). They 
commute with the Hamiltonian as: [h, s^] = kujnSk- Let 
X be the eigenfunction of h with the eigenvalue e then the 
function §kX shifted by k Landau levels is also the eigen- 
function of h with the eigenvalue e + ku>u- These func- 
tions are orthogonal: (x*SkX) = for k ^ 0. Consider 
the matrix of size 2N + 1 reducing the Hamiltonian l|24(l 
in the basis exp(in<p), where n = —N..N. Apart from few 
boundary eigenvalues all middle eigenvalues are equidis- 
tant. The trace of the Hamiltonian 124J1 in this basis is 
obviously zero. The energy shift of the boundary eigen- 
values is limited provided all \ipk\ ~ L Therefore the 
energy shift of middle eigenvalues is bounded by 1/N. 

Let Xj i n ) be the eigenfunction of the Hamiltonian l|24|) 
corresponding to the middle eigenvalue ju>H in the basis 
of functions \n) = exp(zn0). These eigen functions can 
be chosen to be real. Then we rewrite the Ea. l|22|l as 



G R n ,(e) = -i 



ut 

sr^ l \ e ~J UJ H)t 

2^Xj(n)Xj(n )e 2 dt. 

3 



(26) 

Using the orthogonal J2j Xj( n )Xji n ') — &nri and normal- 
ization (XjlXj) — 1 conditions we find that the density 
of state for Eq. lj2UJ) is given by Eq. p0)l . 
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FIG. 1: Two diagrams: a and b, that together give the lon- 
gitudinal conductivity in the limit of vanishing Landau level 
mixing. Dashed line is the short range potential impurity 
line. Full lines are electron propagations with the Landau 
level index being indicated. 



Conductivity in model S. The conductivity ten- 
sor in the model S has been calculated in both cases of 
the largaiS, and smalU^ Landau level mixing. The Drude 
conductivity depends on the transport scattering time 
Ttr- &aa = e 2 riT tr /m, where aa is either xx or yy. Delta 
correlated disorder potential S{q) = 2iru is a special one 
with the quantum scattering time r = r tr . But in gen- 
eral r < r tr - In the absence of magnetic field the usual 
ladder diagrams account for the difference between r and 
T tr - In relevant magnetic fields: a aa <C cr xy , the Drude 
formula reads a aa — e n/wjjTt r . The diagram in Fig. la 
gives the conductivity e 2 n/ui 2 H T whereas the diagram in 
Figl.b gives the transport correction. For ^-correlated 
impurities the diagram Fig. lb is zero whereas for long 
range disorder the two diagrams Fig. 1 cancels each other 
in the leading order. 

Conductivity in model L. We use the Keldysh 
method to find that the longitudinal conductivity is ex- 
plicitly a real function: a xx (£l) = cr(— Q)*. We sepa- 
rate it into the reactive and the inductive parts: Z(ft) = 
R(il 2 ) + ittL(£l 2 ), and we will not consider the inductive 
term here. For the reactive part we find: 



de f(e)-f(e-n) 



,2tt 



n 



where /(e) is the occupation number of states with energy 
e. The effective conductivity of this state reads: 

oo 

a aa (il, e) — 2J Vn/n (ImGw (e - fi) 



lmG n 



,_x(e) + ImG„„,(e)ImG„-i,„'-i(e _ f})), (28) 



where the imaginary part of the Green function is: 
ImG nn /(e) = ,(e) - G# n ,(e))/2. We use the dia- 

grammatic method^ to expand the disorder average in 
the effective conductivity Eq.JSHJ) into a series of Uk- 



a xx (n,e)= f r ReTr(L>(i)n t L>(s)fl 

J J— oo J — oo 



iet + i[e — Q)s 



u(t + s) 2 



E 



u 2 dtds^,(29) 

Z7T 



where U(t) = exp(—iht) is the unitary evolution oper- 
ator of the Hamiltonian l|24|) . The operator II is the 
raising Landau level operator with the matrix elements 
II nin _i = \fn and IF is the Hermitian conjugated opera- 
tor. The integrals with respect to t, s from — oo to give 
advanced Green functions and the integrals from to oo 
give the retarded Green functions. Next we make use of 
the fact that Xj( n ) is the eigenfunction of the Hamilto- 
nian (|24() with the eigenvalue ju>H ■ We rewrite Eq. i|29[l : 



/OO />OG 
/ 

J.K -°° J - 



-<s(l*l+M) 



— oo -J — oo 



■ / \ . / . . s u(t + s) 2 

i[e - juj H )t + i(e-U- ku H )s 

e ^ dtdtyO) 

where the bare dissipation term with < <5 — ► is added 
for the integral convergence. The sum over j and k is re- 
stricted: j — k ± 1. The integral over t + s gives the den- 
sity of states whereas the integral over s gives the Dirac 
delta functions: S(fl ± luh)- Therefore if the magnetic 
phase is neglected then only one of the two Green func- 
tions in the Kubo-Greenwood formula for conductivity: 
a — e 2 uImG£>ImG is broaden. Equivalently the diffusion 
constant in the Einstein formula: er(£l) = e 2 i/(ep)D(£l), 
is that of the ideal electrons: £>(£!) ~ <5(£1 ± ojh)- To 
account for the magnetic phase in the model L would 
probably require a non-perturbative approach. 
In the same way we find the Hall conductivity: 



a xy (fl) 



de a x v(e,n)-a x v(e,-n) 

7T /_„ 2lT 



(31) 



where the Hall conductivity of the electron state e is: 

oo 

a xy {e,9)= ^ v 7 ^ (R*Gn-i,n'-i(e + fi) 

nn' = 1 

ImG w (e) - ReG„„' (e + n)ImG n _ 1 , n /_ 1 (e)), (32) 

where the real part of the Green function: ReG„„< (e) = 
(G^„, (e) + G^ n , (e))/2. We expand the impurity average 
in Eq. l|32ll into a series over coefficients Uk and the result 
is partially summed up as: 

a xv (e,n)= f f f lmTi(u(t)&U(s)fl\sgn{s) 

J J — oo J — OO 

%et + t(e + il)8 2^— d 2^ h 

2 k 2 dtds ±lt m 

2tt k 1 

In terms of the eigenfunction of the operator U(t) (|23|) : 
a^(c,n) = ]T ((Xj&XkHXkflXj) - cc) r dtds 

' J— OO 

u(t + s) 2 



sgn(s) e 



i(e — ju>ii)t + i(e + il — ku>H)s — 



(34) 
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Now we make use of the matrix element: {Xj^Xk) — 
Vjfij,k+ii where 5j.k is the Kroneker symbol. Taking 
the integrals over t and then s we find: a xy (e,n) = 
irp(e)/(u>H + 0)j where the density of states: 



P(e) = 



^exp 



(6- JU) H f 

2u 



(35) 



The Hall conductivity is given by the Ea. (|31|l : 



a xy (ft) 



H 



n 2 



| /<«)*) 



(36) 



The integral over e in Eq.(j2H|) gives the total number of 



At Q = wc find a xy — 



filled Landau levels: N — 2ttI h ti. 

cr H = ecn/H. This result explicitly does not depend on 
the electron distribution function /(e). 



IV. TWO LENGTHS DISORDER MODEL 

Scattering rate off impurities is proportional to the 
scattering cross-section and the velocity of incoming par- 
ticle. The velocity of 2D electron in magnetic field is 
zero. It is disorder that creates the dispersion e(p) - Lan- 
dau level broadening - and imparts random velocity to 
electrons: v = de(p)/dp. In the previous Section the long 
range disorder was shown to be ineffective in broaden- 
ing the Landau level locally. Thus we add a weak short 
range disorder to overcome this problem. No microwave 
excitation is considered in this Section. 

The origin of short range disorder in GaAs heterostruc- 
tures may lie in the barrier layer Gai- x Al x As. Because 
Al atoms are distributed randomly the local energy bar- 
rier fluctuates as ~ 1/ V ' Ni rn where AT,- m is the typical 
number of Al atoms in the volume of typical electron 
wave function 100 A The short range disorder potential 
w(r) is assumed Gaussian and to have the correlation 
function: S w = (w(q)w(—q)) — 1/2ttu(6f)t w — l/mr w . 
Using the vertex Ea. (|10|l . the definition Ea. (|I4l) and the 
orthogonal and normalization properties for Laguerre 
polynomes we find the 'short range' impurity line in N- 
th Landau level: w — Wjstn — S w /2wl H = loh/2ttt w . We 
assume that w <C u, but w could be comparable to Ut r - 

In the limit Ih -C 2d <§; R c the conductivity of the two 
lengths disorder model could be found. We neglect the 
diagrams with crossed short-short and short-long impu- 
rity lines whereas the magnetic phase for crossed long- 
long impurity lines is assumed zero. The average Green 
function is determined by the equation: 



G R {e) = -i 



J exp (- ^t 2 + iet - iZ R (e)t^ dt, (37) 



where £ (e) = wG R (e) is the self-energy of electron in 
local long-range potential before averaging. Solution of 
Ea. (|37J) satisfies the normalization condition for any w: 
C 00 ImG A (e)de = n 



The longitudinal conductivity is given by diagrams in 
Fig.l where the impurity line is necessarily the short 
range line. Remember that the unitary rotation U(t) 
Eq. (|23|) eliminates for long range potential all transitions 
between Landau levels such as in Fig.l. The diagram 
Fig. lb with the short range line is zero. Using unitary 
rotation Eq. (|23J) we simplify the two current vertices and 
the Green functions in the diagram Fig. la and find the 
longitudinal conductivity: 

2N 

cJaa = e 2 ^w (ImG R (e F )) (ImG R (e F )) . (38) 

This equation coincides with that in the phenomenolog- 
ical force-force correlation approaches, but we derived 
the limits of applicability of Ea. (|38l) . In the quantum 
limit of non-overlapping Landau levels loh S> \fu at 
higher magnetic fields: 



Pa 



4 - 



(e F - Nujh? 



e z nT,„ 2ttu 



(39) 



whereas in the classical limit u>n <C \pu in weaker mag- 
netic fields we find the Drude formula: 



Paa = m/e 2 nT u 



(40) 



Corrections to Eas. H39l40|l could arise due to a possible 
additional contribution to the long range disorder origi- 
nating from very long range technological imperfections 
at qR c -C 1. Accordingly we distinct the technologi- 
cal and the donor disorders: u = Uoo + U2d- In the case 
Uca ^> U2d there will be a second crossover at loh ~ y'S 
from the classical Drude regime Ea. (|40(l at weaker mag- 
netic fields to a classical positive magneto-resistivity: 



Pa 



UJ H 



e 2 n,T w y/2nu 2 d 



(41) 



at higher magnetic fields: y/u2d -C ojh *C ^/u^. Because 
Ih 2d <C i? c , \/u2d ~ "vH and, therefore, p xx ~ \fll - 
such positive magneto-resistivity behaviour is often seen 
on experimental traces. 

Drifting states. For long range potential 2d ^> R c 
states of electron transform into drifting states along the 
energy levels and during a short time locally the long 
range potential could be approximated by uniform elec- 
tric field crossed with magnetic field. Crossing of 'long 
and short range' impurity lines now has magnetic phase 
close to zero and one has to account for multiple cross- 
ings to build up the magnetic phase. Formally this can 
be described using the multi-field Hubbard -Stratanovich 
decomposition of the action averaged over the long range 
potential that includes the impurity line weights and the 
magnetic phase factors: 

S[4>] - ~\ J dtdt' S(q)V nin2 (-q)V n3n M 

V4, Wm (-9)Vn 2 (*) V4 3 W)Vm (*') (42) 
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n 



m 



n 



FIG. 2: The self-energy. Dashed line is the short range po- 
tential impurity line. Full lines are electron propagation. 



where the time t is the Schwinger proper time in Eq. i|2U|) 
and operators Vm act on the Landau level degeneracy 
index. For simplicity scattering that does not change the 
Landau level number N is considered. The momentum 
integral in the action is performed symbolically: 



dtdt' ^{t)^{t') 

-1/2 



WMt), (43) 



where A is the area of system, : ... : stands for nor- 
mal ordering and operators a, b act at times t,t' corre- 
spondingly. Next we use the condition qR c -C 1 and 
assumption of short time electron propagation, with the 
expectation values of a and being much smaller then 
2d I In. Therefore we expand the action l|43|) : 



«1M = 4»S(| 

i,j=0 v 



2i+2j 



dijWdntf) (44) 



where for short times the operator O mn has the asymp- 
tote: O mn (t) = C nm ipHt) (a t ) m (a)" ip(t), where C nm 
are combinatorial coefficients. Next we make the Hub- 
bard -Stratanovich transformation using fields 4>q, <pi 



S[rf>]=-i / Hir dt 



(45) 



where A — 1 and the effective long range Hamiltonian is: 



H, r = V50t L Q + l 2L 1 (^ a t + fta) + i,. 



(46) 



Note that Eq. (|45|l is an analog of the mean- field (saddle) 
approximation in the usual impurity scattering theory. 

Next we find the Green function averaged over the 
short range disorder. 'Short range' impurity lines do 
not cross due to the vanishing magnetic phase factor 
and the self-energy is exactly one diagram in Fig. 2 
averaged over the long range disorder eventually. In 
the mean-field approximation for the long range disor- 
der Ea. i|46[l the magnetic phase factors of the two ver- 
tices in Fig. 2 are entangled with the electron propaga- 
tion: (0| exp(i(Q*a + (3a' l ')/v / 2) exp(itHi r ) exp(-i(Q*a + 
Qa))/y/2) exp(—itHi r )\0) . Calculating this factor we find 
the self-energy: 



E(t) = J wV 2 N (Q) e 4d 



tQ* + rtQ)d 2 Q 

2tt ' { ' 



Using the vertex in terms of Laguerre polynome Eq. i|10|) 
we find: £(i) = w (L N (x 2 /2)) exp(— x 2 /2), where x = 
y/u\4>i\t/2d. In high Landau levels N ^> 1 important is 
the evolution time such that y/ut -C 1 and \/2Nx 1. 
Thus we use the large argument asymptote for Laguerre 
polynome: £(£) = wJq(\/2Nx), and Bessel function (ne- 
glecting fast oscillations): 

For retarded Green function t > and further solution 
of the Dyson equation proceeds exactly as in Ref. 8 : 



1 



G R (e, 



— iA log 



8 



Ml) 



(49) 



In the limit to < uR 2 c /(2d) 2 we find ImG R (t, cj)i) = 
exp (— At \og(£/A)). The semi-circle density of states 
transforms into the Lorentz density of states. The av- 
erage over long range potential Green function is: 



(50) 



This result could be derived in alternative way. In high 
Landau levels scattering for electron drifting in the long 
range random potential u(r) reads: 



1 

27 



= to / ImG(e - Qyd y u)\V NN (Q)\ 



. d 2 Q 
(2tt) 2 ' 



(51) 



where we used Landau gauge with y axis being along the 
long range electric field. Evaluating Ea. l|51(l we find: 



1 w (R c 

— = =5 lOg 

r Tr 2 R c \du\ \2d 



(52) 



Diffusion between drifting states. We use the Ein- 
stein relationship and calculate the diffusion constant 
from the diagrams in Fig.l. Averaging over very long 
range potential qR c 1 is done using the mean field ap- 
proximation Eq. I|45|1 . We parameterize Green functions 
by the Schwinger proper times t and s and write down 
the conductivity for example for the diagram in Fig. la 
(after the unitary rotation Eq. l|23|l that eliminates the 
Landau level mixing - m = n): 



VFia.Xa — 

7T 



>2tt tt 



d 2 Q 



V 2 N {Q) (0| Im 



dte -itH lr [a\a] e 



Im r ds e -^Hir[a\a] 



(2tt)2 



V2 



10X53) 



The integral with respect to t + s is convergent in very 
short times ~ 1/ and gives the density of states that is 



8 



finally omitted from the diffusion constant. Then we set 
s = —t and calculate the matrix element in the space of 
boson operators a 1 * x a. It is identical as for the self-energy 
above. We do the same for the diagram in Fig. lb and 
sum up both. Finally we find the anisotropic diffusion 
constant with two components along and perpendicular 
to the local long range electric field (</>i): 



where the rate 1/r is given by Eq. (|52|l with the local field 
\dU p \ ~ U p re^/L 2 . After scattering with the probability 
Eq.(J5BJ electron will hop onto adjacent trajectory that 
diverges from the old trajectory after passing the saddle 
point by distance L p . Therefore we estimate D mac ~ 
CL 2 w/U p > D mic , where C = \og(L p /R c ) log( v / 7rni? c ). 



Di, t = Nw 



dt 



P(t) 
S(t) 



ImG R (t,<j>i) (54) 



where P{t) = J§(£t) - Jf(£t) and S(t) = J (£t)J 2 (£t) + 
J 2 (£t). Using Eg. (|49l> we evaluate the 'diffusion' Green 
function as: lmG R (t) = exp (— At log(£ /A)). Diffusion 
constants for a given local long range field (f>\ are: 



D, 



Di 



l) 



2Nw (\ 



tt£ 



2Nw 
ir£ 

-E(k) 



1 N A 2 k 
k E ^ - 4^ 
A 2 k 



K(k) 
28 



(55) 



where K(k) and E{k) are Jacobi elliptic functions with 
modulus k = 2£/ v / A 2 L + 4£ 2 and A L = Alog(£/A). In 
the limit of weak short range disorder: w -C uR 2 /(2d) 2 , 
the diffusion constant becomes isotropic: 



D = 



2dR r 



TTWU 



(56) 



Thus we derived that for very long range disorder the 
diffusion constant is inversely proportional to the average 
absolute value amplitude of the long range potential. 

Escape from classical localization. Each short 
range scattering event transfers electron orbit by typical 
distance R c resulting in the effective microscopic diffu- 
sion constant D m i C ~ L p R c w/U p Eq. l|56f) . where U p and 
L p 3> R c are the typical amplitude and length of the 
long range potential. The average over electrons macro- 
scopic diffusion constant D mac is estimated as the hop- 
ping length L p squared from one 'localization lake' to ad- 
jacent 'localization lake' over the time it takes for electron 
to diffuse to the percolating level t ~ L 2 /D m i C . Thus, 
Dmac ~ L> m i C . This conclusion is wrong. Let us prove 
that electrons with drifting orbit close to the percolating 
level have much larger diffusion constant and dominate 
the average macroscopic diffusion constant. In the vicin- 
ity of percolating saddle point long range potential is 
approximated as 



U p {r) = U p (x 2 - y 2 )/2L 2 p . 



(57) 



We use parameterization x — r cosh(</>) and y — r sinh(</>) 
with Jacobian: dxdy — rdrd(f>. For <f> 3> 1 we average the 
scattering rate Ea. <|52|) over those trajectories that are 
closer than distance R c to the percolating line: 



1 



9{R C - re - *) rdrdcj), 



(58) 



V. DIFFUSION AND RYZHII CONDUCTIVITY 

Quantum states of an electron driven by uniform mi- 
crowave ac-field are explicitly time depend. There ex- 
ists a unitary transformation into the oscillating reference 
frame where wave functions of electron become time in- 
dependent whereas the disorder potential becomes time 
dependent. The correlation function of short range dis- 
order potential reads: 



S w (t,t';q) = S w (q) e 



iq(R(t) - R(t') 



(59) 



where R(t) = (X(t),Y(t)) is the classical elliptical tra- 
jectory of charged particle in crossed constant magnetic 
field and microwave electric ac-field E(t) = E cos(u>t) 
with frequency lo and polarized linearly along x axis: 



X(t) = 



£Ih u>h cos(wi) 



*2N w 1 - u 



H 



Y(t) 



£Ih u 2 ^ sin(ojt) 



(60) 



where the amplitude of ac-field is conveniently expressed 
in the energy terms: £ = eER c . For long range disorder 
we neglect the time-dependent phase in the correlation 
function Eq. (|59|l due to the small wavevector q. For the 
short range disorder this phase is essential and we trans- 
form the correlation function further: 



S w (t,t';q) = S w (q) ^ J m (qTl)e 



imuj(t — t ) 



(61) 



m— — oo 



where J m (x) is the Bessel function and the size of the 
elliptical trajectory is: 



£Ih u H \Ju 2 +u 2 h 
V2N uj{uj 2 -w 2 h ) 



(62) 



First we take average over the short range disorder for a 
given realization of the long-range potential. This calcu- 
lation follows identically that of Ref. 8 : 



f(e)- f(e + muj-n) 



Ae 2 N 



•E 



de J P„ 
2^ 1 S„ 



ImG(e)ImG(e + mu - fi) (63) 



where the axis x is along and the axis y is perpendicular 
to the ac-electric field and 



Sm= / J m (qK) [Jo(qR c )J2(qRc) + Jl{qRc)} 



d 2 q 
(2n) 2 
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FIG. 3: Diagram for the average current in the left vertex. 
It is a linear response to the inhomogeneous electron density 
with the density gradient being indicated by the arrow. 



Jim) [4{qR c ) - Ji(qR c 



d 2 q 



(64) 



Next we average Eq. (|63|I over the long-range disorder po- 
tential using the method of the previous section. In the 
limit Ir -C 2d <C R c each Green function in Ea. (|FJ5|l is 
averaged separately. The dc-conductivity £1 = reads: 



® XX i @yy — 



Ae 2 N 



Pn 

S, r 



ImG(e + muj) 



(f(e + mw) - /(e))|lmG(e) - J ImG(e) \ (,5 ) 



where the Green functions are averaged. Apart from 
the conductivity anisotropy and the transport correc- 
tions in our P m and S m Ea. (|65|l coincides with the dc- 
conductivity found in Ref£. 

Diffusion constant in driven state. We use the 
Keldysh method to calculate the current in the non- 
equilibrium state of electron system driven by microwaves 
with long-range gradient of electron density. The current 
is averaged during the time much shorter then the long 
range relaxation time. The local current density reads: 



j(f,t)=jG+-(ft;ft). 



(66) 



Diffusion has an origin in the short range disorder scatter- 
ing therefore we expand the Green function in Eq. (|66|l in 
series over 'short range' impurity lines. One of the Green 
functions in each diagram of such expansion represents a 
variation of the occupation number distribution Sf(e): 



8G + - = SfG A - G R Sf, 



(67) 



whereas all Green functions to the left are retarded and 
all Green functions to the right are advanced. They are 
summed up into full Green functions and only one 'short 
range' impurity line remains. In weakly non-equilibrium 
Fermi liquid the variation of the occupation number dis- 
tribution is localized in the energy domain in the vicinity 
of the Fermi level. Therefore we use ansatz: 



Sf(e) 



1 df 
v{£f) de 



n^exp^gr*) 



(68) 



It has a quasi-classical property: f Sf(p) d 2 p/(2ir) 2 = 
n^exp^qr). In the perpendicular magnetic field we ex- 
pand exp{iqr) = 1 + iml 2 H q x j + 0(q 2 ) in series of small 



q, neglecting the magnetic shift operator that does not 
contribute to the current. A typical diagram that ap- 
pears in the result is shown in Fig. 3. It is proportional 
to qhq' and therefore gives an additive to the linear re- 
sponse diffusion equation: j = D\7n. The dashed line 
in the microwave driven state is given by Ea. (|61|l . The 
analytical expression for the sum of all such diagrams 
(n = m) gives the diffusion constant: 



Dxx ^ Dyy 



AN 



de / P m (w) 



2n \ S m {uj) 



|lmG(e)ImG(e + mw). (69) 

Because ImG(e) is positively defined function - it rep- 
resents the density of states - the diffusion constant is 
positive D aa > 0. It coincides with the second term in 
the second line of the total conductivity Eq. (|63() that we 
call the diffusion conductivity <j d . By definition they are 
related by the Einstein formula: a^ a = e 2 v(eF)D aa - But 
the Einstein relationship between the total conductiv- 
ity and the diffusion constant is violated in the systems 
driven by microwaves due to the Ryzhii conductivity (the 
first term in the second line of the total conductivity 
Eq.®): 



R 4e 2 N ^ 

xx±Vyy = —^ W 2^ 



R 



de 
2^ 



(/(e) - fie + muj)) ImG(e)— ImG(e + mu) 

de 



(70) 



Therefore we have proved that the total dc-conductivity 
Ea. (|rj3T) is the sum of the diffusion conductivity and the 
Ryzhii conductivity: a(W) — (Jd{W) + (Jr{W), for any 
intensity of the microwave ac-field. 

How zero resistance states emerge. When the 
power of the microwave excitation provided kuju < uj < 
{k + \/2)ll>h is increasing the electron density distribu- 
tion (originally due to the donor potential) become more 
and more even because the negative Ryzhii current moves 
electrons from the places of lower electrostatic energy 
(higher density) to places of higher electrostatic energy 
(lower density). The correlation function of screened 
donor potential is changing from the completely screened 
Eq.(J7| to the unscreened Eq.@- And the smaller is the 
wavevector the more prominent in Eq.JSJ) a harmonic 
of the donor potential becomes. The average harmonic 
magnitude is q dependent Eq.QJ: 



u(q, W) 



2Tre 2 



\q\ + 2ne 2 v(e F ){a D {W) - a R {W))/a D {W) 

. (71) 

We conclude that the harmonic with q = will diverge 
first as the microwave intensity and a R (W) both grow. 
The magnitude of the long range potential grows to in- 
finity at the critical microwave power W c : <r R (W c ) = 
a D (W c ). At powers higher than the critical W > W c 
the harmonics in Eq.(7IJ that belong to the shell \q c \ = 
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a(W — W c ) are divergent and therefore the local electric 
field will become infinite as well as the amplitude. The 
inelastic electron relaxation due to phonons will proba- 
bly prevent the gradient of the long range potential to 
grow to infinity. q c could be considered as an order pa- 
rameter of an emergent large gradient long range poten- 
tial. Due to the obvious frustration in the electron sys- 
tem, that should handle large amplitude harmonics with 
ever changing wavevectors as the power grows, some non- 
equilibrium glass-like pattern of long range potential will 
probably freeze out eventually. In this situation the diffu- 
sion is dominated by electrons drifting in very long range 
(qR c -C 1) random potential with very large local elec- 
tric fields. As discussed in the previous Section the diffu- 
sion constant will be very small - inversely proportional 
to the magnitude of the electric field on the percolating 
contour. Therefore the conductivity is almost zero and a 
zero-resistance state will emerge. 

On the other hand the Hall conductivity in 2DEG 
with almost even electron density and in the long range 
potential is determined by the fast drift (without clo- 
sures) by electrons that are close to the percolating level 
in the direction perpendicular to the applied Hall elec- 
tric field. This mechanism gives the Hall conductivity 
&xy = enc/H, with n being the electron density on the 
percolating energy level. 



VI. CONCLUSION 

For 2DEG in magnetic field under the microwave ex- 
citation we have found the effect of the Ryzhii current 
on the Coulomb screening of long range external poten- 
tial. The screening factor is divergent at some critical mi- 
crowave intensity. At this and higher powers very large 
amplitude harmonics of very long range potential devel- 
ops. We have found the diffusion constant for quasipar- 
ticles in such large amplitude long range potential due to 
the scattering on weak short range disorder. The larger 
is the long range potential the smaller is the longitudi- 
nal conductivity - therefore the states with almost zero 
resistance emerge. 
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